Self-compensation in manganese-doped ferromagnetic semiconductors 
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We present a theory of interstitial Mn in Mn-doped ferromagnetic semiconductors. Using density- 
functional theory, we show that under the non-equilibrium conditions of growth, interstitial Mn 
is easily formed near the surface by a simple low-energy adsorption pathway. In GaAs, isolated 
interstitial Mn is an electron donor, each compensating two substitutional Mn acceptors. Within an 
impurity-band model, partial compensation promotes ferromagnetic order below the metal-insulator 
transition, with the highest Curie temperature occurring for 0.5 holes per substitutional Mn. 



Ferromagnetism in dilute magnetic semiconductors is 
generally believed to be mediated by carriers — electrons 
or holes — originating from the magnetic dopants them- 
selves. For example, an isolated Mn impurity in GaAs 
can substitute for Ga and contribute one hole, which is 
weakly bound to its acceptor core [Q. GaAs samples 
with Mn dopant concentrations in the range 1-2% are 
ferromagnetic insulators, while samples in the range 3- 
6% are ferromagnetic metals In the metallic phase 
the nominal hole concentration, p, is in principle equal to 
the number of Mn atoms per unit volume. Measured hole 
concentrations are much smaller, by factors ranging from 
-3 for MnGaAs @ to '-lO or more for MnGe §. In most 
theories of ferromagnetism in dilute magnetic semicon- 
ductors, reduced hole concentrations suppress the Curie 
temperature |@, |[ ^. The reverse scenario — raising the 
Curie temperature by increasing the hole concentration — 
is therefore of great current interest. For current theo- 
retical reviews see Refs. ||, |6|. 

Recent experiments show a strong correlation between 
Curie temperature, carrier concentration, and the frac- 
tion of Mn found at interstitial sites Q]. In this paper 
we address several questions not yet settled by experi- 
ment. (1) By what mechanism are interstitials formed, 
given that their calculated formation energies are consid- 
erably higher than substitutionals? (2) What determines 
the relative abundance of interstitials and substitution- 
als? (3) Under what (doping) conditions do interstitials 
act as compensators? (4) What role does compensation 
play in the ferromagnetism? 

To answer these questions, we use density- functional 
theory (DFT) to establish the following: (i) During the 
MnGaAs growth, Mn adatoms follow a very simple low- 
energy pathway to directly form interstitial Mn near the 
surface, (ii) The deposition of additional As converts 
some of these interstitials to substitutional sites, (iii) The 
remaining interstitial Mn atoms act as donors, each com- 
pensating two substitutional acceptors. Finally, we show 
that within an impurity-band model, ferromagnetism be- 
low the metal-insulator transition is most favorable — in 
the sense of the highest Curie temperature — for 0.5 holes 
per substitutional Mn. 

For MnGaAs grown by MBE, recent channeling 
Rutherford backscattering experiments show that as 



much as ~15% of the total Mn may be interstitial 0| . An 
open theoretical question is how interstitial Mn might be 
formed, under what conditions, and in what concentra- 
tion. For a system in thermodynamic equilibrium, the 
concentration of each impurity species is determined by 
its formation energy. Our calculations (described below) 
show that the equilibrium concentration of interstitial 
Mn is negligible. But thermodynamic equilibrium can 
only be achieved if local metastable configurations can be 
readily overcome. We find that the energy barriers sep- 
arating interstitial and substitutional configurations do 
not satisfy this criterion. Hence, equilibrium thermody- 
namics is not a reliable guide for studying the formation 
of interstitial Mn. 

To investigate the incorporation of Mn under non- 
equilibrium conditions, one must identify specific reac- 
tion pathways and calculate their energy barriers. At 
low Mn concentrations, the growth of MnGaAs at low 
temperatures is governed by the potential-energy sur- 
face for individual Mn adatoms adsorbing and diffusing 
on the GaAs surface. We confine our attention to the 
GaAs(OOl) surface, the standard orientation for grow- 
ing MnGaAs. Under As-rich conditions, the GaAs(OOl) 
surface can have several different reconstructions, all of 
which contain surface As dimers as a common building 
block. Thus we studied adsorption on a chemically rea- 
sonable model surface consisting of five layers of bulk 
GaAs plus a dimerized As top layer; the bottom layer 
was passivated. We calculated total energies and forces 
using pseudopotentials and the generalized-gradient ap- 
proximation |8[ ||] . To identify the lowest energy adsorp- 
tion site we placed the Mn adatom at various surface 
sites in a (2x2) supercell, then relaxed the height of the 
adatom and the positions of all atoms in the top five lay- 
ers of the surface. By far the most stable adsorption site 
was the As-dimer bridge site, i.e., equidistant from the 
two As atoms in a dimer. 

To identify specific adsorption pathways, we calculated 
the total energy on a grid of Mn positions in a plane nor- 
mal to the surface and containing the As dimer. The re- 
sulting potential-energy surface and several possible ad- 
sorption pathways are shown in Fig. 1. Beginning from a 
height of several A, Mn adatoms are attracted to the sur- 
face and at low kinetic energies will be "funneled" into 
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FIG. 1: (color) Potential-energy surface for Mn adsorption 
on GaAs(OOl), plotted in a plane normal to the surface and 
containing the As surface dimer. The minimum energy ad- 
sorption site is the subsurface interstitial site labeled i; the 
corresponding surface geometry is shown (light gray for As, 
dark gray for Ga, yellow for Mn). Typical adsorption path- 
ways funnel Mn adatoms to this interstitial site (heavy curves) 
or to a cave site, c (light curves). Inset: Binding energy of 
a Mn adatom centered on the As-dimer; for comparison, re- 
sults are also shown for a Ga adatom. When additional As is 
deposited, the metastable Mn site, s, becomes more favorable 
and leads to partial incorporation of substitutional Mn. 



one of two adsorption channels. The primary channel 
leads to the As-dimer bridge site. As shown in the inset 
to Fig. 1, there are two locally stable adsorption posi- 
tions above the As dimer: a shallow metastable site 2 A 
above the surface and a stable site 2 A below the surface 
and more favorable by 0.8 eV. These two minima are sep- 
arated by a small barrier of 0.2 eV, associated with the 
opening of the As dimer, which is easily overcome at nor- 
mal growth temperatures. (In comparison, Ga adatoms 
also have two locally stable adsorption sites, in agreement 
with Ref. However, for Ga the interstitial position 

is very unfavorable, and the stable site near the surface 
layer ultimately leads to completely substitutional incor- 
poration [0.) Isolated Mn adatoms that are steered 
into this primary channel will generally reach their equi- 
librium position below the As dimer. This position corre- 
sponds to a bulk interstitial site with four neighboring As 
atoms. With the Mn atom at this position, the As dimer 
above it remains intact with a slightly strained bond. A 
secondary adsorption channel leads, with no energy bar- 
rier, to the cave site between adjacent dimer rows; upon 
continued growth of GaAs this site will also correspond 
to a bulk interstitial site. Thus we have identified a set 
of very low energy pathways which initially steer isolated 
Mn adatoms to interstitial sites. 

Although the energetics of Mn on clean GaAs leads to 
incorporation only at interstitial sites, the subsequent de- 
position of additional GaAs can change this situation and 
lead to partial substitutional incorporation. To illustrate 
this, we computed the change in Mn binding energy re- 
sulting from an additional half monolayer of Ga or As on 
the surface. In the presence of either a Ga or As adlayer, 



the surface As-dimer bond breaks and the metastable 
Mn site above it then corresponds to a substitutional 
site in the zincblende lattice. The adlayer completes the 
fourfold coordination of this site and hence may change 
its stability relative to the interstitial site. The change 
depends on the type of adlayer: for a Ga adlayer, the 
Mn interstitial site remains 0.8 eV more favorable than 
the substitutional site, but for an As adlayer this differ- 
ence is reduced to just 0.1 eV. This reduction reflects the 
crossover from the case of Mn adatoms on clean GaAs 
(which favors interstitials) to the thermodynamic limit 
of isolated Mn in bulk GaAs (which favors substitution- 
als), and suggests that during the growth Mn will be 
incorporated at both interstitial and substitutional sites. 

The electrical activity of a Mn impurity is determined 
by its formation energy as a function of its charge state. 
We used DFT to calculate the formation energies of sub- 
stitutional and interstitial Mn using supercells containing 
54 atoms, with convergence checks using 128 atoms. We 
fixed the chemical potentials to correspond to the As- 
rich, Mn-rich conditions normally used in growth. Fig. 2 
shows the resulting formation energies, as a function of 
the Fermi level, for charge states that are stable within 
the GaAs bandgap. Substitutional Mn is an acceptor, 
with stable 1— and neutral charge states. The theo- 
retical acceptor ionization energy is 100 meV, in good 
agreement with the experimental value of 113 meV [0. 
Interstitial Mn is a deep donor, with stable 2+, 1+, and 
neutral charge states; its donor levels have not been mea- 
sured experimentally. For p-type material, corresponding 
to the Fermi level near the valence band edge, each inter- 
stitial will compensate two substitutionals, in agreement 
with previous calculations ||l2|, |l^. The resulting hole 
concentration is p = {x — 2y){4/cr), where x and y are 
the number of substitutionals and interstitials per Ga 
site, and a is the the GaAs lattice constant. 

The consequences of partial compensation for ferro- 
magnetism are different above and below the metal- 
insulator transition. In the metallic phase the Mn- 
induced impurity states broaden and merge with the 
GaAs valence band. Most current theories therefore con- 
sider the carriers as free holes in the valence band. In this 
picture the Curie temperature is proportional to p^^^, 
and hence compensation can only reduce the Curie tem- 
perature. But the role of compensation is less obvious 
for Mn concentrations between 1 and 2%, where ferro- 
magnetism persists experimentally even though trans- 
port data clearly show the samples to be insulating 
The nature of this insulating state is not currently un- 
derstood. If the Coulomb U is larger than the impurity 
band width, then this state is a Mott-Hubbard insula- 
tor and an impurity-band model with localized carriers 
is an appropriate starting point p[ . Below we show 
that within such a model, the compensation affects the 
magnetic interactions in a surprising way: the maximum 
Curie temperature is obtained not for zero compensation. 



3 




0.0 0.5 1.0 1.5 

Fermi level (eV) 



FIG. 2: Formation energies of isolated Mn impurities in GaAs, 
for different charge states as a function of the Fermi level. 
Inset: Transition levels (crossing points of the formation en- 
ergies) between stable charge states of substitutional and in- 
terstitial Mn impurities. 

but rather for a partially compensated system with 0.5 
holes per substitutional Mn. 

We derive this result starting from a Zener double- 
exchange Hamiltonian for the motion of holes in a narrow 
impurity band [p^ , 

= ^ ejol^ajo- + ^ tijal^Uja - ^ • Si. (1) 

zcr ij(7 i 

The first term describes fiuctuations in the on-site impu- 
rity levels caused by the random distribution of charged 
{A~) and neutral {A°) Mn acceptors; the second term de- 
scribes the hopping of holes between Mn spins; and the 
last term describes the exchange coupling between holes 
and Mn spins. The energies Ci are distributed within the 
disorder-broadened band width, W; and aia are cre- 
ation and annihilation operators for holes with spin a at 
site i; Uj are hopping integrals; Jh is the exchange cou- 
pling constant; Si and Si are the operators for the hole 
and Mn (5* = 5/2) spins, respectively. 

It is well established that the coupling between the hole 
and the Mn spin is antiferromagnetic, Jh < 0, and that 
5| Jh| is in the range 100-250 meV |@, This is compa- 
rable to, or even larger than, the acceptor ionization en- 
ergy. Hence the hole at site i is always antiparallcl to the 
Mn spin, and the Hamiltonian can be reduced to a spin- 
less form with the hopping renormalized to tij cos(0y /2), 
where dij is the angle between Si and Sj . 

To second order in the hopping, Hd^ can be mapped 
to an effective Heisenberg Hamiltonian describing the in- 
teraction between Mn spins. 



where rii is the number of holes at site i (thus, ni=0 and 
1 correspond to A~ and A" acceptors, respectively). Two 
features of Hes are noteworthy. First, the effective cou- 
pling between spins falls off exponentially with their sep- 
aration, Tij, because the hopping integrals have the form 
Uj = to('^Ji) exp(-ry7ao), where to{r) oc r/ao for hydro- 
genic wavefunctions and ag is the effective Bohr radius 
of the Mn acceptor Second, because of the factor 

{rii — Hj), interactions occur only between Mn acceptors 
in different charge states, that is, between A'^ and A~ . 
This implies that ferromagnetism is absent in the limiting 
cases of no compensation and complete compensation. 

This result can be easily understood within a simple 
toy model. Consider the Hamiltonian -STde applied to two 
Mn acceptors that share either zero, one, or two holes. 
In the case of zero holes (both A~ acceptors), there is 
no interaction at all. For one hole (A° and A~), the in- 
teraction is ferromagnetic, because the hole can lower its 
kinetic energy via hopping only for aligned spins. For two 
holes (both A'^), the interaction is always antiferromag- 
netic and scales as tfj/Jn. In the real system, aligning 
the spins in the absence of compensation drives the Fermi 
level into the gap between spin-up and spin-down states, 
making ferromagnetism unfavorable. 

To find the Curie temperature of the ferromagnetic 
phase transition, we first note that mean-field theory is 
not well justified because of the exponential dependence 
of the integrals tij and because of the Mn positional disor- 
der |l^ . Instead we represent the system as a percola- 
tion network of randomly distributed sites, each occupied 
by either A'^ or A~ , with interactions only between sites 
with different charges. In the standard solution to this 
two-color percolation problem the critical percola- 
tion radius, Rc, depends on the solution of the one-color 
problem and on the maximum eigenvalue of the connec- 
tivity matrix for the two-color problem. In our case this 
leads to 
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where Be — 2.7 is the average number of bonds within 
the critical radius for the one-color problem and q 
is the number of holes per substitutional Mn. In gen- 
eral, q must be found from the electroneutrality condi- 
tion and the energy levels in Fig. 2, but at reasonable 
temperatures the approximation q = 1 — 2y/a; is ex- 
tremely accurate for y < x/2. Note that due to the factor 
■\/ q{l ~ q), the critical radius diverges for both y = and 
y — x/2, consistent with our earlier conclusion that fer- 
romagnetism vanishes at these limits. 

The Curie temperature is given within percolation the- 
ory by 

kBT,^-Bj^ to{Rcfexp{~2Rjao), (4) 
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thereby also ferromagnetic coupling — can take place not 
only between A'^ and A~ but also between A'^ and 
acceptors. The percolation radius is still given by Eq. ||, 
but q now must be replaced by the fraction of neutral 
acceptors, which must be calculated from the density of 
states. For a simple triangular density of states, the re- 
sulting Curie temperature is shown in Fig. 3(b). One sig- 
nificant new feature arises: ferromagnetism now persists 
even without compensation. Otherwise, the main fea- 
tures of the Mott-Hubbard case are obtained here as well, 
and the maximum Curie temperature is again reached for 
q = 0.5. 

This work was supported by the Office of Naval Re- 
search, the DARPA SpinS program, and NSF Grant 
DMR-0071823. We thank C.S. Hellberg and I.I. Mazin 
for helpful discussions. 
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FIG. 3: Dependence of Curie temperature on hole concen- 
tration for 1.5% Mn, in two different insulating regimes, (a) 
Mott-Hubbard insulator, U > W . (b) Anderson localization, 
with the arbitrary choice U — 0.6511^. Upper and lower Hub- 
bard bands are shown. 



where ks is the Boltzmann constant. The angle brackets 
denote the most probable value among atoms within the 
critical radius. We estimate this value using a hydrogenic 
impurity wavefunction with effective Bohr radius oq — 
7.8 A ImI, and take the energy separation near the Fermi 
level to be the Coulomb gap jl^, given by |e,; — ej \ — 
e'^/nrij, where e is the electron charge and k, — 10.66 is 



the dielectric constant of GaAs 
temperature simply as 



pj. This gives the Curie 
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— ) exp(-2i?,/ao). (5) 
ao 



We show in Fig. 3(a) the dependence of the Curie temper- 
ature on the number of holes per substitutional Mn and, 
equivalently for the case of pure self-compensation, on 
the fraction of interstitial Mn. The maximum Curie tem- 
perature is reached for q — 0.5 holes per substitutional 
Mn, corresponding to an interstitial fraction y = x/A. 

The above theory ignores the orbital degeneracy of 
the impurity band. This simplification is justified if the 
Coulomb energy, [/, of two holes in different orbital states 
is larger than both the hopping integral and the impurity 
band width, W. For the case when U is smaller than W, 
the resulting overlap of the two Hubbard bands will lead 
to a finite density of states at the Fermi level, even with- 
out compensation. For Mn concentrations between 1 and 
2% the materials are nonetheless insulating, suggesting 
that the states are localized by an Anderson transition 
induced by strong disorder. In this case, hopping — and 
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